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ments over which the terms of the heat balance 
are computed. The control mechanisms which 
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PROGRAM DESCRIPTION 


A. IDENTIFICATION 
Program Name 

Bioengineer Contact 

Progranmer 

Date of Issue 


Thermoregulatory Model of Stolwijk (Transient 
Version) 

D. J. Grounds 

V. J. Marks 

June 30, 1974 


B. GENERAL DESCRIPTION 


This user's guide and program description updates and replaces 
TIR 741-MED-3011 . The model is designed to simulate the transient 
response of the human thermoregulatory system to thermal inputs . The 
model consists of 4l compartments over which the terms of the heat 
balance are computed. The control mechanisms which are identified arc 
sweating, shivering, vasoconstriction and vasodilation. 


C. USAGE AND RESTRICTIONS 

Machine and Compiler Required 
Peripheral Equipment Required 


Univae 1110 TSS and FORTRAN V 

Time Sharing Terminal 
Tektronix 4010 Graphics Terminal 


D. PARTICULAR DESCRIPTION 


A detailed description of the basic model is given in references 1-5. 
The current modifications include an improved shivering mechanism and 
the introduction of a graphics output option. 

The shivering mechanism which appeared in earlier versions was des- 
cribed by the following equation. 

QSHIV = CCHIL*C0LD (l) + SCHIL*C0LDS + PCHIL*C0LD (l) *C0LDS . 

But, since CCHIL = 0.0 and SCHIL = 0.0, the expanded equation becomes 

QSHIV = 12.22 (Tcore set - Tcore) head* ( 1 K(l)*(Tset(l) - T(l) 

where 1 equals all skin compartments and K(l's) are constants based on 
the weighted mass of each skin compartment. This equation was based on 
the widely accepted data of Benzinger and Kitzinger(6) . However, the 
data are presented with average skin temperature as an independent 


2 


variable. This does not allow verification of the model equation 
because of the introduction of set point temperatures for each skin 
compartment. The improved modeling of this mechanism is complicated 
by several factors. First, although the role of hypothalmic tempera- 
ture signal and skin temperature signal has each been established as 
necessary for shivering, the functional relationship for the combined 
effect of skin receptor signal with hypothalmic signal is unknown for 
conditions other than uniform skin temperature. Further, the heat that 
is generated by chemical reactions as a response to cold, chemical 
thermogenesis is initiated by the same control signals and is therefore 
extremely difficult to separate from shivering thermogenesis. Shivering 
can also be initiated through the higher centers as a general sympathetic 
response. Also, there is considerable variation in the shivering res- 
ponses between male and female subjects. These factors have contributed 
to the lack of an accepted control system model of the shivering mech- 
anism. 

In order to make the best use of available data without making assump- 
tions which seriously limit the accuracy of the model, the approach 
chosen was to derive an empirical relationship from the experimental 
data. This approach was also taken by Riggs (7), but acceptable corres- 
pondence to the data could not be obtained using his equations. Although 
this approach seems to fulfill the requirements of the model, it should 
be reviewed whenever a more complete understanding of the mechanisms of 
this process is gained. 

By using a combination of linear and parabolic regression schemes based 
on least squares criterion, a set of equations was found to approximate 
the curves of reference 6, page 649. These equations are: 


XTC * TC - (.1 ((37.0 - TC) 1.7) 2 ) / 10) 

RM = 22221. -6lU .2 (XTC) + TS (-1933-2 + 53-66 (XTC) ) 
+ TS 2 (46.45 - 1.289 (XTC) ) 

Q3 HIV = RMS-QRASAL >0.0 


where 

Q3 HIV = Rate of heat produced 
QBASAL = Basal Metabolic Rate 
RM = Metabolic Rate 
RMX = Metabolic Rate 
TC = Intercranial Temp. 

TS = Average Skir Temp. 


by shivering 
BTU/HR 
CAL/SEC 
BTU/HR 
°C 
°C 


BTU/HR 


( 1 ) 

( 2 ) 

(3) 


The fitted curves are shown in Figure 1 with corresponding experimental 
data. 
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DESCRIPTION OF INPUT 


4 


The user must specify if the graphics output capability is required. 

If graphics is required, then a list of variables to plot is solicited. 
Upon completion of the list a name list input is requested. This 
request is answered as shown in the example of Appendix 2 by supplying 
the following values: 

Col. 2 


RM 

Total Metabolic Rate 

BTU/HR 

QBS 

Basal Metabolic Rate 

BTU/HR 

VEFF 

Efficiency of useful work 

i P 

AC 

Surface Area for Convection 

FT % 

AR 

Surface Area for Radiation 

FT 

TCAB 

Temperature of the Cabin 

°F 

TW 

Temperature of the Walls 

°F 

TDSWC 

Dew Point Temp, in Cabin 

°F 

VCAB 

Cabin Free Stream Velocity 

FT^/8EC 

VEFF 

Ventilation Efficiency 

i 

PCAB 

Atmospheric Pressure in Cabin 

PSIA 

0 

Gravity Normal to Earth 


CLOV 

Clothing Thickness/Conductivity 

FT 2 HR °F/BTU 

EUG 

Emissivity of Undergarment 


CPG 

Specific Heat of Gas 

BTU/lb°F 

DT 

Integration Step Size ( .05) 

MIN 

PR INTI 

Print Interval 

MIN 

SETT I 

Max Time /Run 

MIN 

MCASES 

Number of Imposed Case 


MCASES = 0 

Initial Conditions (Not plotted) 


MCASES > 0 

Imposed Case Using Final Values 

for Starting Point 


To complete the namelist, input a $ END is entered beginning in column 2. 
F. DESCRIPTION OF OUTPUT 


See Appendix for example output . The graphics output option which has 
been added to the program is demonstrated in the example run shown. 

It includes some conversational input to specify the variables to plot 
and their ranges. There are 55 variables which can be used with the 
graphics option: 


1-43 


Temperatures of the T Array 
(See TIR 741-MED-3011) 

°F 

44 

TUGAV 

Average Temperature of Undergarment 

°F 

45 

SQUG 

Sum of Sensible Heat Transfer Rate 

BTU/HR 

46 

QEVAP 

Heat Transfer Rate of Evaporation from Surface 

BTU/HR 

47 

TOTL 

Total Latent Heat Transfer Rate 

BTU/HR 

48 

TDEWC 

Dew Point Temp, in Cabin 

BTU/HR 

49 

WORK 

The Metabolic Work Performed 

°F 

50 

QSHIV 

Heat Produced by Shivering 

BTU/HR 

51 

STORAT 

The Rate of Heat Accumulation 

BTU/HR 
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52 

QSTOR 

Total Heat Accumulated 

53 

TCAB 

Temperature of the Cabin 

54 

VPDEW 

Vapor Pressure at Cabin Dew Point 

55 

U 

Useful Work Performed 


BTU/HR 

BTU 

°F 

BTU/HR 


G. INTERNAL CHECKS AND EXITS 


The model will proceed until a steady-state condition is met (Storat - 
Oldstor/Storat) <..01 or until maximum time is exceeded whenever 
MCASES - 0. The model is then initialized with the steady-state values 
and thermal transients are then simulated for MCASES >0. 


H. INDEPENDENT SUBROUTINES 

The Independent Subroutines are MANT, SHRT, VPT, CQSH, TCF, GVAR, and 
the Tektronix Plot Package. The Subroutines are listed in Appendix 1. 

I. SYSTEM SUBROUTINES 

No special system routines are required. 

J. COMPLETION DATE 


June 25, 197 1 * 


I 
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c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CNO 

* T Ol»tjK METABOLIC Haw TRANSIENT model 
COMMO»J/NSH |V/TAVG,TC,RmX 

COMMON T 1 9 3 I , T UG A V . S Q'JG • QE V AF » ToTl .TDERC.FOijK , QSh I V , 
t ST0RAT»QST0R,TCA8.VP0EN»U»TUGI 10) ,ACE( 10) iArCIIOI t 
6 PCftA,RH,QLCG,C(91),TSET(93l,TR,£UG,CL0.C"G,vEFF,G. 

(. VCAR ,0T JpE iwd»S»i 

COMMON /O I AO/QCONv < 90 » , QC0N0t9D> ,BFt9t>) .QHET(rO) »TEST I 91 t ,rArh[ri 1 - 
DIMENSION SA VTUG I 10) t S A V T (1)1 i PC A I 10) » PLT ( 55 ) , N 1 (8) 

FOuIVALENCE «PLTU»*TI|»» 

DATA PCA/.Q7,.3602,.0670S,.Q67nS,.ISB7 l .ISG7,.02S,,025,2».Q393/ 
DATA iTSCTl I » • |«l |93)/9S.6.97.6 ,9?.0 , 66 . 6 , 96 . 9 , 98 . 3 iff . 9 , 99 . 9 , 94 . 1 , 
. «S*1,9n. I, 93 . 7,96.1, ?5 M ,9N. 1,93.7, 97. 6.94*5, 

• ~ 9*. I ,’9.5,97.6,96,5,95. | ,99. 5,95*9, 95.7, 95. 6, 

• 9 5.5,9S.9,9s.7,95.6,95.s,95.8,9S*5,9S.7,9S.5, 

• 9s.«t, 9S, $, 95,7, 95, 5, 96 . 5 , 95,, 9$,/ 

OATA C/6. 67, 1.938,0. *,96, 0.535, 70.9,39.9, 9, 92, 2 . 69,1.3*2 
. »3» 35, 0.699, 0,989, I, 382, 3 . 35, 0«699, 0.9*9, 9 . 9 , 10. 2 

• .1.59,1. 192, 9,9, 10 . 2,1. Sail, |92, 0,1568, 0.0738, 

0*0992,0. 1 89 ,0. 156* ,0.0738 ,0,0992,0.189 ,0.2695 , 

. 0. 0738 , (», ih9,o,297,0. 24^5,0.0736,0.198,0.297 , 

. 9.79/ 

OATA »r»Ra/lHy,0.«2/ 


DEFINITION of BODY SEGMENT TEmpERATuRE SUBSCRIPTS 


T 4 I » 

• 

ME Ay 

COME 

— — — 

T 4 2 ) 

T I 9 ) 

■ 

head 

SK 1 N 


T ( 5 ) 

T ( 7) 

n 

trunk 

FAT 


T 1 8 ) 

T 1 1 0 ) 

9 

9 I GMT 

arm 

muscle 

T ( 1 1 ) 

T 1 1 3 > 

m 

left 

ARm 

C 0 H Y 

T< 19) 

T ( 1 6 ) 

u 

left 

ARP 

SK 1 •« 

T ( I 7 ) 

T< 19) 

9 

pjGmt 

LEG 

FA f 

T4 20 ) 

T ( 22 ) 

9 

left 

leg 

MUSCLE 

T ( 2 3 ) 

T ( 25 ) 

9 

RIGHT 

HAND CURE 

T ( 26 ) 

T ( 28 ) 

9 

RIGHT 

HAND SKIN 

T ( 29 ) 

T 1 31) 

9 

left 

hand 

FAT 

T( 32) 

T ( 39 ) 

9 

right 

FOOT MUSCLE 

T ( 35 ) 

T 1 37 ) 

m 

left 

foot 

CORE 

T ( 3 ft ) 

T « 90 ) 

9 

left 

FOOT 

SK I N 

T ( 9 l ) 

T 1 93 ) 

m 

average MUSCLE 



HEAD 

MUSCLE 

rt Jr- 


MEAD 

FAT 

trunk 

CORE 

T ( 6 ) 

■ 

TRUNK 

MUSC 

TRUNK 

SK 1 N 

T I 9 ) 

■ 

RIGHT 

ARM 

RIGHT 

ARM f AT 

T ( I 2 ) 

■ 

R I GmT 

ARM 

LEFT 

ARM MySC L E 

T< 15) 

• 

left 

ARM F 

RIGHT 

LEG C°«E 

TUB) 

■ 

R 1 GmT 

LEG 

R l GMT 

LEG S*lN 

T 1 2 1 ) 

• 

left 

LEG C 

LEFT 

lfg fat 

T ( 2 9 ) 

■ 

LEFT 

leg s 

R 1 GMT 

HAND Md S cleTI27) 

• 

RIGHT 

hANO 

left 

HANO C 0Re 

T C 30 ) 

■ 

LEFT 

hand 

left 

hand s*In 

T 4 3 3 ) 

• 

R I GMT 

foot 

RIGHT 

FOOT FA T 

T 1 36 ) 

■ 

RIGHT 

FOOT 

LEFT 

FOOT m U5cl€ 

T< 39) 

• 

LEFT 

foot 

central blqoo 

T ( 92 ) 

■ 

AVERAGE SK 


NAmEl 1 ST / INPUT /RM»«b 5 ,tlEFF, AC I AR I TCA8 ,T*, T 0C /C ivCAB,v£FP, PC AB 
0 i<..C;_OV,EUG.CPG,OTtPRlNTI,sEri,MCASES 

COMMON /PL T BUf /NBuF , A0«*F 4 »P| I ,VBUF< 181 ,fl» ,nA4 3» iKSTOPPiSETI iNN 


*R I TE I 6 , 1 0 > 

IT FORMAT ( / • 000 YOU *AMT GRAPHIC INSTEAD CF tABuLAR OUTPUT? (Y/N)»» 


READ 15.15) KYY 
1 5 FORMAT < 1 A 1 ) 

IF (KYY.EQ.KY) KPlT- 1 
IF (rpt.T.GT.O) CALL £v ft R 4 N 4 r N ft^ NW 1 
2 I IF (KPlT.EQ*! . ft N U , NBUF.GT.O) KSTOPP-I 
IF ( ICONO.EQ.OJNKPLT.kPL t 
IF (KSTOPP.LT. I ) GO TO 25 
call PLOTX 
KSTOPP-O 


UKlUiNAl r 
OF POOR QUAl 


N6OF*0 


2 S AR I TE ( 6 » 30 ) 


transient Thermal mooEl 
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JO FORMAT ( /‘OINPJT OATA uS«M6 NAMELIST »|NPUt*'/» 
RE AD ( S* INPUT ) 

Qd ASAL-QBS 

IE ( QhaSAL.LT .0.01 > UBASAl«2«J* 

IE (QBA5AL.6T.MM) g ft a S A L • R M 
*« I TC U .INPUT ) 

IE ( MCASES.E4.0 > ..RITE »6,A0; 

4a FORMAT l l*C* I 7*. IdWlRj T lAt CONDITIONS//) 

I) (MCASES.NE.m ..HITE <»»8C> 

AD FORMAT I iHUi 1 7X » IdHlMPnSEO CONDITIONS//) 

DT I Mr«DT/60. 

SET*«SET|/AC. 

IE(NKPLT.E0» » .ANO.HCASES.GT.OKPlTmI 

SiTT.sETX 

PR I WT«PR INT 1 /AO. 



c 

C INITIALIZATION 

c 



PMN0'*'»PR I NT 

If ( mcasES.NT .0 ) 60 TO | mQ 

IE (ICOND.EU.I .AND. MCASES.EO.OI WRITE«A,6Sm) 
ASM FORMAT ( ' NOT 1 14 1 T I A l l 7E0 AT STEADY STATE*) 

kpl t»q 
ICOND-0 
F 0 UN T ■ 0 
OLCG-O, 

SOUG«0. 

QEVAP-O* 

totl-o. 

QSh I van. 

OSTOR-O. 

SToRAT-O, 


00 100 I • 1 . ** 3 
1 0 (- T l I JaTsETC I ) 
TUr,A VaT ( M2 » 
PO 120 l*I»U) 

J"Ma I 

120 TU6I | ) a T ( J ) 



I Ml DO | AC I a 1 . 1 0 


M • I 

ACE ( I ) aP c A ( 1 ) a AC _ 

ARE ( I ) aPCA ( I ) *AH 
1 AC CONTINUE 

ST I ME»ST I mE*T IME 
TlMEao. 

U«UEFF/ 1 OC.*)RN-hrAsAL) 
WORKaRM-^BASAL-U 

VPDEWbvPT ( TOE PC ) 

CLOa , B8*CC0V 
l*>C CONTINUE 
2uC FORMAT ( 1 2SH 1 T IhE 

|0 EVAP Q LATENT hEAT 

1 - - -*/~ L2SH MIN 

3 R T u / H R 8TU/HR STORAGE 


ORDINAL p. ... 
OP POOR QUALJt 


Temp av te«p temp q 

SHIVER tot;l cabin de* 

HEA& Sp I N , ONOFIT*- SENSIBLE 

rate HEaT temp POINT 


TI.AN5ICMT TmIRMAl M00EL 
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B / 

S 

A / 

I2SH 


CORE 

f GaBMCNT BTu/MR 

V b M 

n B 7" 

BT*'/7*B 

F 

BTOrABe t 

F 

7 

LC«7 


B TU/HR 


Br«//) 


i tf < a » 78* > 

ruKH«T(iH i lx • • ?m i n* • , 1 1 x • *HCi>nco*c avB'xI* tugav blood 

CVAP I^StNS S T 0 s A T QSMlv GS TOR fCAB TDERC*) 

C PA| N L0 0p FOR SHlKTSLfFVF C A 5F 

• • 

221' P T | M aT|M£«6U. 

IF (LC.6T.SH) GO TO 1R0 
LC«LC*l 


IT (PCASFS.EQ.OJGO To S A 7 
NBtJF «\ilUF ♦ I 

IT t N bl)F . G T * I 8 1 ) N»iUF"l*)l 
XtJjF I NRUF ) "P T I M 
00 2J0 1*1 ***< 

Nl 1»NI ( I ) 


nEfii ^ 

<**?£* qoa ^ 


23f, YBUF ( NBUF , I » "PL T ( J I 1 » 

IT (FPlT.GT.O) CALL PI.OT* 


IT (rPLT.GT.cn GO To 2*0 
S 6 7 CONTINUE 

C RHlTE(6.2»TAVG*^fi» .TC.OSHIV 

2 FORMAT(Ih ,FI I »%'£*, »NfT R A T«. • • f t T« S , 

GSX.M-C TEMP** »f I » . b ,5X * »OSH I V«* ,F | | ,5 ) 

C U0R87 I L • l . 3 3 » 8 

C 8 8 7 .tHITE(A»8BB> < 0C0N V ( J L- I ♦ I LL > • I|.L* I »8J 
888 FORMAT ( 8 f9»2> 

I TE ( A • 2** 0 > P T 1 M • T ( I ) , T • N 2 ) i T >JQ A V , T ( R | ) ,Q r VAp,TOTL iSTORAT • 
6 *SM I V i 9ST0H • TC AB » TOEWP 
2 Nn F0RMAT(F8.||9X»IIF9.2) 

26l> OL 05 T R • STONAT 


QE 


C 

C 


CALL smut 

...... — 

T 1 ME«r I ME .QT I HE 

IF ( ABS ( STOHAT ) • GT . 2 . ) 60 TO 280 

IT ( AqS< (OLOSTH-STOHAT I/STORAT ) »lT. .001 » r, 0 T 0 300 
28 U IF ( TJmE.GF.SETT ) GO TO 300 
IT ( MC ASES.EW.*) > GO TO 2*0 
IF ( PRNOw * GT . T I ME ) GO TO 260 

prnot»prnor*print 
GO TO 220 
300 pT I Ms T I M£»*Q * 

P«NO*»T I ME.OT | ME 
- if (MCASES.EG.O)GO M-~S*4— 
nBuF»NBUF. 1 

IF ( MBUF .GT • I 8 I ) NBUF • I B J 
XBUF ( NRUF ) »PT i M 
00 310 1*1 fNM 
H I I ■ N I ( I I 

31*' V&tiff NBUF , | l»PLT ftff H- 

IF (KPLT.GT.O ) CALL plotx 


TRANSIENT ThIRHM. HOO^L 
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If t«Pl.T,GT«OI «.*. TO 3 1 S 
SA* COnTIhoC 

v»H|TC |A,2**0)PT |h,T( I | ,n«lt) tTuCAVtTlRl I » qE V 4 P . TOIL , 5 T OR A T • W5h I V . 
•WSTOR.TCAB.TLf *C 
3| r * K0ONT aK OUN T ♦ 1 

IF I IC0*0*£U* I •ANO.KOUNT.LE.Mf ASCS) GO TO JO,, 

C 5 *vE INITIAL CONDITION RESULTS To START Im^OsED CONQ|T|On CASES 

32U DO J4.3 la* # *l3 

S A y T I I I ■ T I I • 

3RD CONTINUE 




00 3 AO I • 1 t I 0 
SAVTUGI I I • T U G < I I 

360 CONTINUE 

A T UG A v • T OG A 9 
XSquG-SQuG 
XLEVAPaQC VAP 

xtotl-totl 

X STR A T*S TOR A T 
X 0 S H I 9aQ3M | V 
XOSTOPaOSTOR 
I CONO« I 
GO TO 20 
C 

C RESTORE INITIAL CONDITION RESETS FOR NEXT Case 

c 

36*' 00 MOO I ■ 1 * ** 3 _ 

Till aSAVT < I » 

**00 CONTINUE 

00 N 20 !■! • 10 

TUG ( I I a S A V T UG ( | J 
•*20 CONTINUE 

- TU6A VaX TUGA V 


SWUGaXSQUG 
QE VAPaXQC VAP 


TOTL axTOTL 

STORATaXSTRAT 

qShiv-xqshiv 

USTOeaXQSTOR 

GO TO 20 

END 

SUBROUTINE SHRT 

COMMON T I H 3 » . TUGA V .SOUG iCEVAP , TOTL # ToE»C ,RO k K ,q 5 Hl V » 
k STORAT tOSTOR , TCAb . VPDEA »U»TUG I 1 0 ) » ACE ( 10 ) * ArE I I 0 ) » 

6 PCaB,RM*nLCG«C Ml ) ,TSCIli 3 )»TA»CUG«ClO*CpG«yEFF , G t 
k VCAR,DTIMF*GB*«SAL 

C0MM0N/SHRTNN/EMAX<l0l,QRSENl,0RsEN2.0RSEN3.GRSENS*QRSEN* t 
• USE it I I 0 I • QR AO I 1 0 ) 

0 I MENS ION H I I U I 

DATA H/«033».02A,2*.036,2*.03 , » t 2a.0**i2»«03A/ 

TARbta*HAU« 

SOUG A aO • 0 
SOUGRaO.O 


SONaOtO 

TAvSKN»l0. 1 *‘*6*Tlb)»n.Ofl26*Tll2l*0»082A*Tl|6) 4 .0,19N5aTl20)40.I9R*v« 

6T I 2R 1 ) / • 9902 


C • • • i 

C 


TRANSIENT Thermal M 00C L 
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C CALCULATION C F W-RA&I ATEP(GRAO) ANO g - 5 E N 5 I B L C ( g S f N > 

• •••• 

CO *0 !•! I |0 
J«*N# I 

U'GRaTUO I | I • N 6 0 • 

MC«H| | > *ACE » I I »S«RT IPC aB»VCAB > 

IFIG.Lr«0.0)CO TO 10 

HtlaO.n«*ACC(I>*(rCAb*«2«<j«A0S(TUG(n-TCAB>)»*.2S 
IF (HC 1 » G T • MC >HC»HC I 

10 M h • 0 • I 7 1 3E“®*ARE t 1 • •Ll'G* ( T U G R • • 3aTuGRaa2*T*RaTUGRaT * R • • 2 ♦ 

|F ( I ,lT * 2 » 0 H • I • GT . A » GO To 70 
ificlo.lt.o«oi »GU TO 70 

TUC( I »*IhA*Ti.*MC«TCAH*ACH ! l/CLO»TI JI»/lH»*Mc*ACCl I l/CLO) 

GO TO NO 

20 TUG( I )«T ( J) 

NQ fiUGRanRa ( TUG I | I -T t. I 

gUGAanC* ( TU6I I )-TCA*> 

J-ttfG«»»S9UG»^<»UG* 

SQuga«souga*luga 
QSfN ( | | • QUG A 
aK AO I I ) agUGN 
AO CONTINU 



c CALCULATION ot hFSPITORY sensible 

c 



ON SEN l a 0,5*0. ON I 6 a PC AO* I * 1 , 0 / I NB. 3 a I TC ABa<|6R # AN ») aRMaCPGa ( C 0* SbSaT 
.(lM0»n^A*T<2)*0.07<i7*T«3)^0#73B*TIS)*0»2#,l&*T<A)l-TCAB» 

«*<SEN2 a p « I 7 2 • NHSt»'1 
QNSEN 3 ■ 0.US7N • QRSCNl 
QRSENA • 0»&23 • LNSENI 
QNSENS ■ 0.*< 7 * • URSENl 
ON SEN | aO, 77 l agHSfcU I 

C 

c 

SG1JG a SqUG* ♦ Sgua.. ♦ S Q * ♦ c R S E N | a Q R S r N * ‘ Q R s F N 2 ♦ Q R 5 1 N Y ♦ 

.QNSE*'A 

TUGAvao»33l 7 *TUG»2»a0»|ONa(TUGl3»»TUG(N))*0»230l5*lTUG<5)aTUGIA)) 



c 

c C A LCtjLATE MAAlMCm EVAPORATION RATE 

c 

00 60 I ■ I • 1 o 
J"N* 1 

vFTUGaV*»T I TUG (II) 

HEaO.I 2 AaACE‘»l««TCAB*NAO.»aa| .ONayEFf / « Oo • *sOPT I VC AB # PC A 6 > 

IF (G.LE»0*0)G0 TO AS 

HL|al.32aACE(I)a(TCABaNAO.I/PCABa(PCABaGa ( ABsl.OOSapcABaTuG(I)- 
aTC AB I a | .02»<VPT<TuG( I I l-WPDEW) ) )»*,25 
|F(mFI • GT • HE I HE«MEI 
AS I F ( I . L T, ? .OR* I *G-. A) GO TO 70 
IFICLO *L T • .01) W T» 70 

HECLa27»36*ACE( l )a(T(J)aN6T.)«a0,8|/(CL0apCAH) 



r\r> n n n r> n n n 


TUANS I EmT TmIRNAL HOOEl 


lU 


C <ii E A T ,SH I vt* iCONSTRI C T | ON ,0!u AT ION CALCULATION* 
C 


c 


c 

c 

c 

c 

c 


DO 00 

TE$T» 1 ) *T • I >-TS€T t I ) 

6 ARN( I )* 0*0 

COLO* n*o.o 

If * TEST ( | ) I 20 iHO ,60 
in COLD ( l ) ••TEST ( I ) 

H<) 60 TO 80 
60 RARMI 1 ) "TEST < l ) 

80 CONTINUE 


POo » <**ny 


T A vsr«T < S2 I 

ti ♦n;s37»*Armfli*T,nm i«a» h<i2I*o«o*ii*"armm 6> 
♦0»093**ARM(2n>»0,093».‘/A'»M<2‘U*O.Ol I0 7 S*AARM<28)*0»0I 1 07S* 
RARH ( 32)+0«0|99S*MARN( 36)*0»0|99S*ffAAN(R0) 

COLDS b O»082 7 »COLD( ** ) ♦0»S87*C0L0( 8 ) ♦ 0 « 0 ** 1 J * C 0 l 0 ( 12 I *0 » 0 M 1 l • C 0 L D ( 16) 
•OOT3*CDLO<2n>to.OT3*COl.Ot2 , <)*0»nno7S*COLDT28)'^D»OIl075* 
C0 l 0(32)^0.U|/9S*C0L0(36)*0*0199S*C0»-0(‘»0» 
RARMH»0«9iT«.iiARMt6)+0.n9S»*ANM( 1 () I *0 • 09 5 • R A R M ( 191*0*1 90S**ARM( I 8 I ♦ 
0»196S»AARH«22) 

C0L0m« 0»9 I 7»C0L0 ( 6 ) ♦0,o95*cOLn ( 10 > ♦0«P95 *c0Lq II H I *0. 1 90S*COLO 1 I « » ♦ 
0» I 96S*COLD « 22 ) 

5*€ AT-C5w*#A*m ( | » *S iw*w A rms ♦F«W*W ARM < H •****5 
DlLAT«C0 IL»* A AO( I »*SDIL*AAWMS*PDIL*AARm(1)*#aRMS 
STrIc»CCON*COLOI 1 » *SCON*COLDS*PCON*COLD I n*C0L0S 
C * LL ClSH(TIlIi'2SHlv,TAVS«.QrtASAL) 


CALCULATION Of RESPIHATOPr LATENT 


QLAT l**5*0«091 9 • P C A H • 1 <4M, / ( M H • 3 * ( TCAB*9S9,69) )*Rm«(vpT(0,38S* 

•T< I » ♦0*086«T< 2> ♦0.0287*TC 3 I *0. 236»T«SI ♦0.2 a 15«T<6M-0.8»vp0E«O 

•• I ( I 8 .0* I 090*0 > / » 3 2 . U«PCAB ) ) 

QLAT7-0* I 72*NLAT I 
QLAT3«,0579*'3LAT I 
QLAT5*0»R76*WLAT 1 
5LAT6»fj»S23*QLAT 1 
QLAT 1 • 0 • 7 7 l *QL A T 1 
QfT»QLATl>QLAT 2 > 0 LAT 3 ^QLAT 5 ^(aLAT 6 


S < I N DIFFUSION 


00 100 1 ■ 1 t 1 u 

ooif H 1"6*66* ♦ VPT < TU6( | » » - VP0F W 3 ♦ AC€ < 13^ 
1 no CONTINUE 
120 CONTINUE 


QcATENT(WLAT) CALCULATIONS 


c 



TRANSIENT TmERM*L MODEL 
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00 140 ! • 1 *10 

QLATIl)*Q0|FlIlaSK|NS<l»aS*EATa2.aaUT(J»-TS £ ;T(Jn/N,) 

IF ( oc A T ( I )*GT«EhAX( I l I 60 TO |N 0 
GO TO | 6 0 

1*0 QL AT I | ) "(MAX ( I ) 

1*3 CONTINUE \ 

T 0 TL »00 

00 | AO !■! t »0 

QDaQoaQD|F ( I > 

| #0 T 0 Tl«TOTl» 0 U AT ( I ) 

QEVAPaTOTL-QR-OD 

C BlOUO FLOW CALCULATIONS 

••aa«aaa.aaa«aaa»aaa aaaaaa 


/ 


c 

c 

c 


00 200 1*1 • 1 0 
N*H * |-3 

^ BF ( N ) «BFP ( N 1 

QH£T 1 *» »QB< N ) *NHASAL 

QMET I Na 1 ) aQB I Na 1 )*GBASaL»hO*K'M I )a*ORKaCN|LM ( I > aQSH J y 
BF(n*1 )*BFB<N*I >^'JH£T(N»i »-OB(N+| >*qbasal 
GHET(N* 2)*QH(N-*2)*QAASAL 
-*f INa2T»HTOiN*2» 

r OMCf fNaJl VQB 44 L+J ) — ■ — \ 

aF-L4^rr«7BFB (N^J l *SK Uv I | ) •£> IL AT I I ,0*SIL t N€ t I »*5?AM-€* — 

200 CONT fNUE 

TSBEaBFI H >aBF ( B | a*F ( l 2 ) *6F (|4> a*F I 20) *BF ( 2 * UBF (20 >a*F <32 )*bF < 3*> 
. ♦BF(NO) 


CHECK FOR NEGATIVE BLOOD FLOA 


DO 220 I ■ I • NO 

221 IF ( BE < I I ,LT »0.0 l BF < l )«n.O 

C 

C OCONV I I I "CON VEC T 1 ON FROM BLOOD TO EACH NOO^ 

c QCONDI I laCONOUCTlON BETWEEN ADJACENT NODES 

DO 2H0 I • 1 i N 0 

QCONV I ! >«BF ( I) • ( T < <H )-T< I I ) 

QCONDI I I*F ACTOH II I • ( Til I -T ( I ♦ I ) ) 

2 Hi CONTINUE 

a • a 

C TEMPERATURE CALCULAT I ONS — - 

Caaaaaaaaaaaaa # aaaaaaaa«aaaaa*aaaaaaaaaaa*aaaaaaaa 9 aaaaaaaaaaaaaaaa 

C ...... ...... ... ... — ... 

C CALCULATE Temp OF head CORE.TIIIi And TRUNK c°RE»T<5). 

C--- .... ... ..... — .................................. 

Tll)aT(II*DTlNE/C<Ua(QMET(I)-oLATiaQCDNV(l).QCONO(l)-QRSENl) 

TlS)aTl5)aOTlME/Cl5i * l QMET l 5 1 -QLATS* QC0NV<51-QC0NDf&)-QRS£NS) 

— ......... 

C CALCULATE TEMPERATURES OF REMAINING CORES — a«M ( 9a l 3 J * cEG ( 1 7*2 1 J , 

C HANDl25a2?) .AND F00T(33«-37 j 

C ....... .... — 

DO 260 I»’i 37, M 

T I I J «TI I iafl TIME /Cl 1 1 aloMETI I 1 ♦ qCON V (4 ) -QC<>MO<I J-L 

2 61 CONTINUE 


u u u 


transient thermal model 
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c c*tcuL*Tr the temp E iMTimrs or the wusctt m«nnrt*t »**Mrup 

c 1*0. LE6(18*22), HAND (26*30), F 00 T(3 m*38) 

C-- -. 

T(2)aT(2)«DT!ME/C(2)*(QCONO(l)*QMCT(2>-3LAT2*QCONV(2)-QCONO(2)- 

• QRSEN2 1 

T« 6 )«T( 6 )* 0 T!ME/CU)*(( 3 C 0 NI 3 (S)*qMET( 6 )- 3 i.AT 6 *QC 0 NV( 6 )- 0 C 0 ND( 6 )- 

1 QFSEN& ) 

00 2 R 0 l* 10 i 38 ,N 

T<I)»T(|)*DT|ME/C(l)*(OCONO(l-l)*OMET(I»*QCONV(t >-9C0NDt I ) » 

280 CONTINUE 

c — - — — - — ------ 

C CALCULATE TEMPERATURES of THE FAT layer --HEAD(J) |TRUN<(7) ,ARM( I |*15) 

C LfCi< |R*23» iHAN0(27*31 » ,FO0T< 3S^3^) 

........ ............................... 

T(3)»T(3)*OTlME/C(3)*(QCONO(2)*OMET(3)-gLAT3*QCONVI3)-«CONO(3l- 

• 3RSEN3) 

00 300 I »7 i 39 , N 

T ( I) aT ( I) * 0 T 1 ME/C < I ) • ( QCONO ( I - 1 ) ♦QMET ( I > *qCOhV ( I ) -QCONO ( I ) ) 

30 " CONTINUE 


C CALCULATE TEMPERATURES OF THE 3*1** --HE AD ( ** ) , TRuH(C < 8 I • AR M ( I 2 ♦ I 6 ) , 
C LE<j(2U*2N) ,HAND(2fl*32) »F00T(36*NQ) 

C— — ........... .......................................... 

DO 320 I-M.RO.M 
J«I/N 

T ( I) »T ( 1 ) *0T IME/C ( I ) • ( oCOMO I I-| ) *OMET ( I ) -qLAt ( J I *QCONV ( I ) 

• -OSEN ( J » -ORAD ( J ) -OLCG ) 

3 2C 1 CONTINUE 

C-- — ......... ............................................ 

c CALCULATE TEMP of CENTRAL M'.OOD --(Ml) 

C-- ........ — ....... ... ... ..... ...... 

SOCONV ■ 0,0 
DO 3 H 0 I • l • MO 
SQCONV«SOCONv-QCONV ( | ) 

3no CONTINUE 

T(Nn«T<MlJTDTIME/C(*U)*SQCONV 


CALCULATE average SKIN temperature < H2 ) BASED ON PERCENTAGE of 
TOTAL SKIN AREA FOR £A C H S*IN NODE • THAT NOo^S TEMPERATURE 


T»H2»«0*07«T(M)*0.3602*T(8)*O.OA70'iPT(l2j*0»o6705*T(16»PO»l587* 

, T( 20 )*G, 1 & 87 *T( 2 M)* 0 ,U 25 PT< 28 )* 0 « 025 »T( 32 )* 0 , 03 M 3 *T( 36 )* 

• 0 . 03 M 3 PT ( MO ) 

T‘M3)»n»0232S*T(2)*0,5R9*T(6)*0*0527PT(tO)*0,OS27*T(lM)*0.1S92* 
, T(18)*0.1SV2*t(22)*0.0OI1S*t(26)*0.00iIS*T<3O) 4 'O.001I5* 

, T(3M)*0,0ul 1 5 • T ( 38 ) 

TBF*n # 0 
DO 360 1*1 • M 0 
3*0 TSF«?ttF*frM-l-» 

POlSE«S»926*TBF/60,0 
RETURN 
END 

SUBROUTINE C(iSH(TIC.Q5Hlv.TAySK»QBASAL) 
common/nshiv/tavg,tc,rmx 

- CHt TCF-f TC » T 1C » O* 

CALL TCF( TAVG.TAVSK ,0) 


TRANSIENT ThEHHAl Moore 
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TC*TC-(*l-(<<37.0-TC)-|.7>*-2.)/|0.) 

RMR-p 2 2*» .-MR.2-! I*33*?*5 3.A8»M r l F - 

&«TAVG«-2*I9*.9S-| • 289* ( TC ) ) 

C CONVERTS NCTAHOLlC RATE f ROM CAL/SEC TO BTU/NR 
RHX-Rnx-3.6'3.97 
QSH I VaRMX-QBASAL 
IF iqshI v.le»o.o>os..i v-o.r 

RE TU-U 

END 

SUBROUTINE TCF (TC.TF.IFc) 

IF C |FL.NE.O*<»0 TO AS 
TC*(TF-*2* >*&•/*• 

r L turn OW5INAL MQt H 

4 S tf.t € **«/%.*32 OF POOR QUALITY 

RETURN 

END 

FUNCT I ON VPT ( T ) 

C FtjNCTION Tq CALCULATE VAPOR PRESSURE AT TEMp-T 
X - A*47.27-(T*‘*AP. )/| .8 

TEMP- A* I .«/ < T*9*0. J • ( 3,29 9*6, 8 *8£-3-X«-l « I 7 OE 3 > 

6 / I I . ♦ 2 • 1 B8E-3* * ) 

VFT»1207./|0.**TENP 
RE TUP I. 

END 

SUBROUTINE GVARlNl . N A . N N ) 

— DIMENSION N 1 I B ) »NA lh | t «_AI BS ) 


data 

IkA( I 

1 • 1- 1 • 

S ) / ' 

T ( 1 ) 

T( 2) 

T1 3) 

T ( R ) 

T ( S ) '/ 

data 

( K A ( I 

) • I ■* ♦ 

1 0 > / • 

T I * ) 

T ( 7 | 

TIB) 

T<9, 

Tl I0> •/ 

data 

( K A ( 1 

) . I * 1 1 

. 1 s » ✓ • 

T ! 1 I ) 

T ( I 2 ) 

T( 13) 

T< l9> 

T( IS) •/ 

oata 

! K A ( 1 

) . I ■ 1 6 

.20/* 

T ( ) * ) 

T ( I 7 > 

T( 18) 

T * 1 9 ) 

T 1 20 ) * / 

data 

( K A ( I 

1.1-21 

. 25 ) / • 

T ( 2 1 ) 

T ( 22 1 

T ( 73 ) 

T * 29 1 

T12S) */ 

- OATA 

4KA4 l 

1*1*2* 

♦ 30) /• 

T ( 2* ) 

T(27> 

T ( 78 ) 

T ! 2 9 ) 

T ( 30 ) */- 

DATA 

< K A ( I 

1.1-31 

,35)/* 

T ( 3 1 ) 

T ( 32 » 

T ( 33 ) 

T<39l 

T(3S) •/ 

DATA 

(FA! I 

1.1*3* 

» NO 1 / ' 

T ( 3* ) 

T 1 37 ) 

T ( 38 ) 

T < 39 > 

T(90» •/ 

DATA 

( F A I 1 

1 . l-Rl 

» 9 b 1 / • 

T ( 9 1 ) 

T ( 9 2 ) 

T ( 93) 

TUGa v 

SQUG * / 

DATA 

1 F A ( I 

) . I -9* 

.50)/* 

©E V AP 

TOTL 

TDE*C 

WORK 

©SHI V • / 

data 

( K A ( I 

) . 1-bl 

. S5 ) / • 

S T OF A T 

QSTOR 

TC AB 

vPUE* 

u*/ 

RRITf!«*IO 

l- 

— 

- - 

— 

— 



— • 

1 1 FORMAT 1 /'OPLEASE 

ENTER 

ITEMS 

TO BE 

PLPTTEo Bv 

INDEX NO. ( J2) • ) 


DO 3D 1*1.8 
GO TO lS 
J 2 A R I T E I A t l R J 

t“ FORMAT!* ••••• READ ERrOK'I 

--lb -READ l S * 20 .ERK-- .L3J *-U4 1 - - 

20 FORMAT ( 1 2 > 

If ( N I C I ) »LT. 1 ) GO TO 90 

IF ( N I ( I ) .GT.SS) 60 TO 13 
M-N I ( i ) 

N A ( I J ar A ( M > 

AR l TFI-A.25) N1 ( 1 ) ,hU444 

2b FORMAT!' •*•••• .I3.2X.A6) 

30 CONTINUE 
Rf NNa | • l 
RETURN 

end 

eno onsite printout cn june 2*. 1979 at lsisojin 
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** f*00» i ^ 

* to* LIff 


APPENDIX 2 
EXAMPLE RUN 


DO TOU WANT GRAPHIC INSTEAD OF TABULAR OU T PU T< > <T/H) 
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•INPUT 



20 


t 


N 

(A 

• 

O 

10 


o 

o 


v> 


<n 


/» 


s n 

IJ 3 • • . 

3UMMA 
K Q O) « ■ 


U» 

X 

o 



(•U • • • 

3 *- » • n 

« n b <n 


♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦♦ i ♦♦ 


o 

U K 

u 

B O 


OKB 

_J 

UBB 

« 

Q 9 • 




ut 


II H 


II 


«• Wt 

■ uu 


« Ui « u. « 
usauyu 


K tat 
3 ZhU 

O V «0 

J3t*-«UU 

if uyuoo.uic 


« • 


t.M .If 'S9.BB M.il 67.fi 


IICASCS"! 

TCAB*4# 

TU-4A. 

•CM® 

•INPUT 
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tel 

<9 


O 

o 

o 


z> 

« u 
(i a 



x 


nr)«< MMMCiMMMM 


♦ ♦♦♦♦♦♦♦♦♦♦ 
telteJtelteiteltelteitelUltelteJ 



♦ ♦ 
tel tel 



tel tel 




a 

x 

o 

u 


w 

o 


X « 

in « 
«u * 
a »- t 
« i 



o ►- 
u m 
a o 


h 


Ifl te. « tel « U « 

CMtelUKUXOUteJO 
kta««H>-te 39 i 


h tel 

a x«« 

OUtf Hte« 

jaiteiuu 

uuuuoimc 



TYPE SMIFT-OUT (SO) AMS RETURN 
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ft 



x 

a 

o 

fr- 

et 


x 

»- 

at 

« 


</> 


A 









A 






% 








• 



• 






l/l 


• 



• 



• 



• 



• 





u> 



u? 



(A 



(A 



(A 



« 


U. 



u. 



U. 



u. 



U. 





<M 



ft 



CM 



CM 



CM 



at 


% 



% 



% 



* 



% 



Ui 


• 



• 














• 



• 



• 



• 



• 



z 


♦ 



♦ 



♦ 











W. 



u. 



U. 



U. 



U. 





% 



% 



% 



% 



• 



Ui 


X 



X 



X 



X 



X 



Z 


<D 



CD 



CD 



• 



® 



•>« 


• 



« 



% 



« 



% 



►* 


♦ 



♦ 






■* 



♦ 



% 


« 

• 


« 

• 


« 

• 


« 

• 


« 

• 


in 



• 



• 



• 



• 



• 


X 













• 



• 

% 


Ui 



Ui 



Ui 



Ui 


% 

Ui 


• 

► 


-1 

3 

• 

-J 

3 

• 

_l 

3 


_! 

3 

• 

_i 

3 

to 



« o 

fs 

ftt o 

to 

<x 

O 

• 

ft o 

(ICO 

♦ 

>- 


u 

-i 

r u» 

—1 

rs 

Ui 

_J 


Ui _l 

1 

Ui 

_i 

t 

r> 


in 



in 



CO 



in 



to 



a. 

















>- 


>- 

Z 


>“ 

X 


>• 

Z 

• 

>- 

x 


>- 

Z 


=> 



a 

• 



• 


«# 



U 



tf 


o 



+4 

<r 


** 

to 


M 

to 


M 

• 



• 




X 

et 


X 

at 


Z 

♦ 


z 

• 


z 

• 

ui 

• 

















• 


u> 



u 



Ui 



Ui 



Ui 


z 

CM 


© 

• 


o 

• 


o 

• 


o 

• 


o 

• 

a 

a* 


-1 

v4 


_i 

ON 


-J 

CM 


_i 

CM 


-J 

CM 

c 

















at 



A» 



AM 











t* 



V) 



in 



V) 



to 



to 



• 


% 



• 



• 



% 



% 



• 


Z 



z 



X 



X 



X 





% 



% 



« 



• 



% 





>• 



>- 


Si 

y- 


H >- 


at >- 




«■ 



CM 

w 


+4 

w 



•a* 


o 

W 



• 


u> 


♦ 

*- 


x 

►« 


at t- 


t- 

fr- 





w 

o 


lit 

o 


o o 


to o 





-i 


i- 

—i 


Cf 

-j 



_i 


a _» 



>- 


ft. 

► 


(L 

>* 


a. 

>• 

to a. 

>- 


a. >* 


A 



A 



A 



A 



A 



A 
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